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Abstract 

We study the one-loop renormalization of high-energy Lorentz violating four fermion models. We 
derive general formulas and then consider a number of specific models. We study the conditions for 
asymptotic freedom and give a practical method to determine the asymptotic-freedom domain. We also 
point out that in some models the RG flow contains "rational" Zimmermann trajectories that might hide 
new symmetries. 
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1 Introduction 



Although Lorentz symmetry is one of the most precise symmetries in nature [T], the possibiHty 
that it might be violated at high energies or very large distances is still open and has been 
extensively investigated. A Lorentz symmetry violation at high energies allows us, among the 
other things, to renormalize vertices that are otherwise non-renormalizable. This result is achieved 
using a modified power counting criterion, which weights space and time differently [2]. Modified 
dispersion relations improve the large- momentum behavior of propagators in such a way that, in 
the common perturbative framework, the theory remains unitary, local, polynomial and causal. 

Using this knowledge, it is possible to formulate a Standard Model extension [3ll3] that is CPT 
invariant, but violates Lorentz symmetry at high energies, and contains two scalar-two fermion 
vertices, as well as four fermion vertices, at the fundamental level. The inclusion of CPT violating 
terms is also possible. Four fermion vertices are important for a variety of reasons. On the one 
hand, they can be used to explain proton decay. On the other hand, they can trigger a Nambu- 
Jona-Lasinio mechanism and give masses to fermions and gauge fields even if the elementary Higgs 
boson is suppressed [1]. In its simplest version, the scalar less Lorentz violating Standard Model 
schematically reads 




where the quadratic terms are 

G 

>Ckin/ = E -jrP' + bi^'P) x'l (1.2) 

a,b=ll=l ^ J 

and the vertices are denoted symbolically, namely without listing all possible field differentiations 
and index contractions. In our notation hats are used to denote time, bars to denote space. The 
weight of time is equal to —1, while the weight of each space coordinate is —1/3. The gauge 
couplings g have weight 1/3. The weighted dimension of space-time is 2, so the Lagrangian 
contains only terms of weights ^ 2. Moreover, Xi = = (z^L'^l)' = Q\ — (^L' X3 = ^°k-, 
Xt = uj^ and xt = d% i^"" = (fg, f/^, J^t), = (e, ^,r), u'^ = {u,c,t) and d°- = {d,s,b). The sum 
is over the gauge groups SU{3)c, SU{2)l and U{1)y- Finally, T = —D'^/K\, where A/, is 
the scale of Lorentz violation, and are polynomials of degree 2. 

The models of O 1^ are anomaly-free, because gauge anomalies cancel out exactly as in the 
Standard Model |3]. The "boundary conditions" such that Lorentz invariance is recovered at low 
energies are that h]°'^ tend to d""^ and tend to 1 (one such condition can be trivially fulfilled 
normalizing the space coordinates x). 
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An important consequence of the high-energy Lorentz violation is that all gauge interactions 
are super-renormalizable, therefore asymptotically free. Moreover, since fermions have weight 
1/2, the four fermion interactions are strictly renormalizable. At energies much larger than 
vectors decouple and the model (jl.ip reduces to a four fermion model in two weighted dimensions, 

a,b=ll=l ^ J ^^L 

plus free fields. The purpose of this paper is to study this type of model, its one-loop beta 
functions and the conditions for asymptotic freedom. We stress that if the high-energy model 
(|1.3p is asymptotically free, then the full Standard Model extension (II. ip is, as well as its other 
versions of ref.s [3lH]. 

We work out a method to determine the domain of asymptotic freedom in quantum field 
theories with more couplings and apply it to some of our models. Our approach is to study 
the asymptotic expansion of the running couplings around the free fixed point. The domain of 
asymptotic freedom Paf is determined by the arbitrary constants contained in the expansion. 
The dimension of Paf is equal the number of positive eigenvalues (including multiplicities) of a 
certain matrix A/", which depends only on the one-loop coefficients of the beta functions. 

Finally, we point out the presence of special RG trajectories that might hide new symmetries. 
Indeed, if we apply Zimmermann's "reduction of couplings" [5] to our beta functions we find that 
some solutions of the RG equations exhibit features that normally appear only in the presence of 
hidden symmetries. 

The paper is organized as follows. In section 2 we classify the four fermion vertices and 
present the most general CPT- and rotation invariant four fermion model. In section 3 we work 
out general formulas for its one- loop beta functions. In section 4 we study some explicit examples. 
In section 5 we formulate our method to determine the domain of asymptotic freedom. In section 
6 we recall Zimmermann's reduction of couplings, and explain why some of our models might 
possess hidden symmetries. Section 7 contains our conclusions. 

2 Four fermion model 

Using charge conjugation we can use only left-handed fermions, which we collect into a vector 
iy We orthonormalize the kinetic terms iji'y^dfj,ik and simultaneously diagonalize the quadratic 
terms ijij'^df^d'^ik by means of a unitary transformation. 

The most general GPT invariant, rotation invariant four fermion vertices are of the form 
L*L*LL and LLLL, plus its conjugate L*L*L*L*. We do not have to include vertices L*LLL. 
Indeed, in the four-component notation they can be constructed only with an odd number of 
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space-time indices, so they violate CPT. Precisely, a term L*LLL plus its Hermitian conjugate 
reads 

a{LiUL2){IJiVU) + a*(Z27°f/7°Li)(Z47V7°L;^), (2.1) 

where U can be 7" or 7*, V can be 1, a^^ or a*-', and a is a constant. However, the combination 
dHJ is CPT odd. 

On the other hand, all combinations of terms LLLL, L*L*L*L* and L*L*LL are CPT even. 
Explicitly, we have the structures 

(LiyL^)(IfyL4), {LiUL2){UUU), 

for L*L*LL, and 

(LfVL^XLlVL^), 

for LLLL, plus their Hermitian conjugates. 

Using Fierz identities, we can show that every LLLL vertex is Lorentz invariant and of the 
form 

(LfL2)(LfL4). (2.2) 
Similarly, all four fermion vertices of type L*L*LL have the form 

(4L2)(4l4). (2.3) 

In general, this structure is Lorentz violating, but the combination 

(l1L2)(4l4) + (l1L4)(4L2) = {LlLl){l4Li) = i(Li7^L2)(Z37^i4), (2.4) 

is Lorentz invariant, which can be easily proved using the Fierz identity (^^ap^^^yS ~ '^^aS^^yii — 

At this point it is convenient to switch to the two-component spinor notation. Because of 
(j2.2p and (|2.3p . the most general four fermion vertices are constructed with the contractions 

a,(3, . . . being spinor indices and denoting transposition. The most general CPT- and rotation 
invariant interaction Lagrangian reads 

^4/ = ^ E [(4eM^m)gijkm + {eje£,)i£le£m)hjkm + {iUt^MstJ f*jk^] , (2.6) 
^ ijkm 

where the couplings / and g satisfy the symmetry properties 

9ijkm 9kmij 9jimki fijkm fkmij fjikm- 



Because of (|2.4p . the most general Lorentz invariant interaction Lagrangian is (|2.6p if the 
couphngs gijkm satisfy the additional symmetry property 



9ijkm — 9imkj ■ 

Finally, the most general high-energy Lorentz violating four fermion model that we are going 
to study reads 

^ = E 4^ + ■ ^||) + ^4/, (2.7) 

where a are the Pauli matrices. This model is renormalizable. The quadratic terms of 
(jl.2p have been switched off, because they are not important at high energies. Nevertheless, the 
correlation functions at generic external momenta have no infrared divergences in two weighted 
dimensions. Therefore, we can study the model (|2.7p in itself. 

Assigning the axial charges +1 to ij and —1 to £|, the ^-terms are axially symmetric, while 
the /-terms explicitly violate the axial symmetry. Suppressing the /-terms we obtain a restricted 
model 

Cr = Yl 4^ + ^^-^ • + ^ E i4mAU9^Jkm, (2.8) 

j ^ ^ ijkm 

that is still renormalizable. This restriction is interesting for merely theoretical purposes. How- 
ever, for phenomenological applications we include the full set (j2.4p of interactions. Integrating 
the Higgs field out in the Standard Model produces /-terms, which generate g(-terms by renormal- 
ization. Moreover, our Lorentz-violating theories can include four fermion vertices that describe 
proton decay, and such terms are of both /- and (7- types [Gj . 

It should be emphasized that a high-energy Lorentz violation does not imply that the proton 
must decay. In this respect, we have two classes of renormalizable models. The models of the 
first class are described by i?-invariant Lagrangians. Then, i?-violating vertices are not generated 
back as counterterms by renormalizatiorf^, so this choice is consistent. The models of the second 
class contain i?-violating four fermion vertices at the classical level. Consistency with existing 
experimental bounds on proton decay imply that the energy scale of Lorentz violation Kl must 
then be greater than or equal to 10^^'^^GeV (see for example [6]). At present, there is no reason 
to expect that A/^ is much smaller than this value. However, if neutrino masses have the Lorentz- 
violating origin suggested by one of us in |3,j, namely they are explained by the renormalizable 
dimension-5 vertex 



^The B-violations due to the B + L anomaly are non-perturbative, and do not affect the renormahzation 
structure of the theory. 
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where H is the Higgs field, then could be around lO^^GeV, which is still compatible with the 
bound coming from proton decajl. 

3 One-loop renormalization 

In this section we study the one- loop renormalization in the most general model (|2.7p . It is 
convenient to use the background field method. Replace £ with £ + -0, where ip denotes the 
background field. Then expand C m. ip and keep only the quadratic part in I. The result can be 
written as 

— X ^ ^ f^rn^n ) (Qn^np + Hnp) I 1 ' 

^ np V > / 

where Qn collects the kinetic terms and Hnp is a matrix quadratic in -0. In momentum space, 



Thus, the one- loop contribution to the effective action reads 

Ti = -^tr ln[Q„(5„p + Hnp] = ^tr[Q'^HQ'^H] + constant + finite. (3.1) 

Observe that the tadpole —iti[Q^^ H]/2 vanishes, since the propagator Q^^ is odd in p. For this 
reason there is no wave-function renormalization at one loop. 

We use a dimensional regularization where only the dimensions of space are continued to 
complex values 3 — £2 • There is no reason to continue also the time dimension, since the integrals 
over p = converge. The divergent parts of just two integrals are necessary to evaluate expression 
()3.ip . namely 

f dpi d^-^^p pI _ Aj . 

J (27r)4 DnD, 47r2e2(|6„| + |6J) ^ 



and 



/ 



dp4 d3-=2^ 1 (p2)3 Ai 

— + finite, 



(27r)4 Al DnDp A7T^e2\bnbp\{K\ + \h 

where 



D.=pl + bj^^ + 6' 



(in Euclidean space). The fictitious mass 6 is introduced to avoid IR problems at vanishing 
external momenta and can be set to zero after the evaluation. The calculations are performed 



difference of one or two orders of magnitude can always be due to the dimensionless couplings that multiply 
the vertices. 
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using Feynman parameters to integrate over p^. This isolates the pole of the p- integral, therefore 
the divergent part. 

Using these formulas we easily find 



^ ^ tnp \ ^^[HnpHp^^ tr 



(47r)2e2 

^ ' np 



a'* I a" i ^ 



We can convert this formula from the dimensional regularization to a conventional cut-off A on 
the p-integral replacing l/£2 with InA+constant. Now, since A has weight 1/3 and, by definition, 
the high-energy dynamical scale ^ has weight 1, matching weights and dimensions we find the 
identification 

1 1 A3 
m ■ 



£2 3 ImK\ ' 

The relation between the one-loop bare and renormalized Lagrangians is 

The beta functions are found equating the /x-derivative of Cb to zero and performing a number 
of straightforward manipulations. We find 

Pg^ ~ 2g^2 ^ ^ ^"P [6(1 ^np)9ijnp{,9pnkm 9pmkn) "^^npOinpjdkpnm 
np 

inkpQpjnm ~\~ ^^fjpmnfinkp) ■^np){9inpm9kpnj 9inkp9pmnj ^^fmpjnfinkp) 

+48 (1 + Snp){fjmnpfiknp ~ fmpjnfiknp ~ fjmnpfkpin)] {i^j)*^[k,m) ' (3-2) 



P f — 1 Q 2 ^ > ^"P [2(^ ~^ ^np) 9nipk fmpjn ~l" 3(1 + Snp)9nipj fkmnp ^Snp9nipkfj\ 



-|^g^2 / > ''"P L^V*^ "npjynvpkj mpjn ~r "J^-^ ~r ■^npjynipj J kmnp ~r '^^npHmpkJ ]pmn\(j^ 
np 



and = 0, where s„p = bnbp/\bnhp\ and t„p = 1/(|6„| + \bp\). The expressions on the right-hand 
sides of (j3.2p have to be symmetrized as follows: 

\-^ijkin\{i,j)ir^(k,m) ~ 1^ {-^ijkm ~l" -^kmij) i 
[^ij'fcm] ' = — {Xijj^fn + -Xjikm ~l~ -^ijmk ~l~ -^jimk") • 



4 Explicit examples 

In this section we consider some particular cases in detail. A separate section is devoted to the 
conditions for asymptotic freedom. 
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U{Nl) model A simple model is the C/(A'^]^)-symmetric model of left-handed fermions L, 

^ L ' L L 

Here we can set 6=1 rescaling the space coordinates, the fields and the couplings. We have a 
model ([221) with 

9ijkm — dl^ij^km ~l~ 92^im^kj j fijkm — 0, 

SO the one-loop beta functions read 

Pi = ^(Nl-2), j3^ = JL-{2Nl-1). (4.1) 
The solutions of the RG equations read 

ffi(i) = ffi(O) + ^^f^ {g2{t) - 52(0)) , 52(t) - 



2Nl - 1 ^ " ' ^ ' 1 - {2Nl - l)52(0)t/(367r2) ' 

where t = — \i].{\x\^) and x is some scale. 
In the ultraviolet limit {t — )• oo) we have 

51 (i) ~ 51(0) - TTITT 752(0) - — — fi(2(i) ' 



INl-I'''' {2Nl-1)H' ' {2NL-l)f 

In particular, 

iVi - 2 , , 

51(00) ~ 51 (0) - 2jy^ _ -^ g2(0). 

If 5i(oo) 7^ the UV fixed point is interacting, if 51(00) = it is free. However, an interacting 
fixed point is not guaranteed to survive beyond the one-loop approximation. Thus, we must 
restrict to the subspace with 51(00) = 0. 

It is easy to prove that the model with 52 = is renormalizable, for example introducing an 
auxiliary field x of weight 1 and writing the Lagrangian in the form 

C'^ = L^H (d + ba- V + xiL^'L') - ^x'- (4.2) 

Such a model has vanishing one-loop beta function. However, it is unlikely to be finite, since no 
symmetry appears to forbid higher-order corrections. 

We have just pointed out a general feature of our four fermion models: there exist combinations 
of couplings that have zero one-loop beta functions. Nevertheless, we are unable to use this 
observation to prove the existence of interacting fixed points. That would require more knowledge 
about higher-order corrections. Using only one-loop results the best we can do is to project onto a 
suitable subspace of the space of couplings, using a method inspired by Zimmermann's "reduction 
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of couplings" [5], and study the conditions for asymptotic freedom in that subspace. We first 
prove that there exists an analytic solution of the RG equations of the form 

oo 

91 = ciff2 + ^ Ciff2 (4.3) 

i=2 

to all orders in the perturbative expansion. Consistence with the RG equations gives 

dffi ^ /?! 

which in turn uniquely determines all coefficients Cj's. We find 

^ Nl-2 
2Nl-1' 

plus recurrence relations of the form 

where Pi{cj^i) are well-defined polynomials depending only on the coefficients Cj's with j < i. 
The Zimmermann solution (j4.3p restricts the two-parameter space to a curve. There our theory 
has a unique coupling, §2, and its one-loop beta function is still given by the second formula of 
(j4.ip . The condition of asymptotic freedom is thus 

92 < 0. 

"Electroweak" model Now we consider a four fermion model containing one family of the 
electroweak model. We have the left-handed doublet L"" = [vl^^l] and the right-handed electron 
cr. The high-energy four fermion Lagrangian reads 

C^W = L^^i (5 + bLa ■ + (^d - bna • cr + -l(LtV^j)(eJjL») 

Because of hypercharge conservation /-terms are not allowed, so this model is of restricted type 
(|2.8p . Define P' = (z^L,eL,e^), where ^ denotes the charge conjugate. Then we have 

gijkm = giL^ij^km + g2LSimSjk + gR^iS^jS^kS^mS 

+ {^ — 9LR){SijSk3^m3 + ^i3^j3^km) + HSimSk3^j3 + 5j3<^m3<Jjfc) > 

where 6ij = 6ij — Si^Sj^. Applying (|3.2p we obtain the beta functions 

PIL - ^^on I , P2L - 3— - + 2— - , Pr 



A / A(s-3) \ ^ A / A(2s + 3)\ 
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where s = 6L6ij/(|6L||&ij|) and 

= ElL _L 2^ -u -E^ 4- ^^^^^ 

"""i&lI i^lI i&Li + i?>ijr 

Recall that bi and 6/^ do not run at one loop. 

We see that the beta functions depend only on three couplings, precisely A, g2L and the 
combination u. The conditions for asymptotic freedom are studied in section 5. There exist 
RG-flow trajectories with special properties, possibly hidden symmetries. This topic is discussed 
in section 6. 



U{Nl) X U{Nji) model This is a generalization of the electroweak model, where the left- and 
right-handed fermions are in the fundamental representations of U{Nl) and U{Nf(), respectively. 
The Lagrangian is 

52 \ . _ / ^2 



where o, 6 = 1, . . . iV^ and /, J = 1, . . . Nr. Define t = (L'^, E"^), with i = (a, /). We have a 
restricted model (|2.8p with couplings 



gijkm — giL^ij^km + g2L5im^jk + giR^ij^km + g2R^im^jk 

+ {g2LR — giLR){5ij5km + ^km^ij) + g2LR{SimSjk + 5 jk^imj 1 

where 5 and 6 are the Kronecker tensors of U{Nl) and U{Nji), respectively. The beta functions 
are 

o g2LR ( , g2LR{s-'?>)\ r, fi'2Li?/ 5'2Lr(2s + 3)\ 

where now 

— EIIl _i_ I /V -I- AT ^^-^ _|_ ^^9iLR 
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Dirac fermions So far, we have considered only explicit examples of reduced type ()2.8p . Now 
we consider a model of N Dirac fermions, which involves also / terms and new types of g terms. 
We impose the flavor symmetry U{N) and parity invariance, which allows us to set bi = bji = 1. 

As usual, write the Dirac fermions as (^1,^2*)) where i\ denote the left-handed components 
and £2 are the charge-conjugates of the right-handed components. The action of parity reads 

pe{ = e£i*, pei = -ei\*. 

Moreover, i\ and £2 belong to the fundamental and anti-fundamental U{N) representations, re- 
spectively. 

The Lagrangian (|2.7p becomes 



91 

2Ai 



92 

2AiL 



(4.6) 



L L L L 



+ 



+ 



/2 



{i{^eii){e^^ei^) + {ifet2n{i\^et2 



and all couplings are real. A possible third vertex of /-type, namely 

is not included, because it is not independent of the other two. Indeed, a Fierz rearrangement 
gives the identity 



{fi'emlUli) + (£fV4)(4^e4) + {^I'eiim^eii) = 0. 
We report only the beta functions in the large N limit, which simplify considerably: 



N{9l + 9l) 



367r2 
N 



/3, 



92 



N{9l+9l) 
187r2 - 



/3, 



93 



Ng29i 
187r2 ' 



/3, 



94 



Ng29i 
97r2 



187r2 

N 
'367r2 



N 



{9I + 9596 + 9I + 12 ff - 12/1/2 + 4/1) , 

{9I + 45556 + 56 + 24/f - 24/1/2 + 4/1) , 

N 
18^ 



= [hig^ + 256) - 3/1 (55 + 56)] , Ph 



■/2(56 - 55)- 



Observe that the beta functions do not depend on gi and g^. Moreover, the couplings separate in 
two groups, (7i_4 and 55,6-/1.21 which do not talk to each other. These, however, are only features 
of the large N limit. 
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5 Asymptotic freedom 



In this section we study the conditions for asymptotic freedom in the presence of more than 
one coupling. The idea is to search for solutions of the RG equations as expansions around the 
ultraviolet free fixed point. The domain of asymptotic freedom Paf is then determined by the 
free parameters contained in the solution. A different approach to asymptotic freedom with more 
couplings is due to Zimmermann [7]. An investigation that is in part related to this problem can 
be found in [8]. 

We first illustrate our method in the case of a single coupling a with beta function 



oo 



/3q = d = Q PnOt^ ■ 



n=l 



where the dot denotes the derivative with respect to t = — ln(|x|//), |x| being some scale. If /3i 7^ 
the asymptotic expansion around the ultraviolet limit t — t- 00 reads 



1 ^ 6n(lnt) 



t ^ t 

n=0 

where 6„ are polynomials of degree n in Int. Inserting this expansion into the RG equation we 
get 60 = — 1 and the recursion relations 

(n - l)hn -h'^ = 6n< (5.2) 

for n > 0, where 6n< is a linear combination of monomials Y[i with ^ • riiki ^ n and depends 
only on the coefficients bm with < m < n. Consider first n = 1, and observe that (5i< contains no 
logarithms, so bi is a polynomial of degree 1. The coefficient of Int in bi is uniquely determined, 
while 61(0) = b remains arbitrary. For n > 1 the relations (15. 2p can be solved recursively: 

-(5n< +bn) = -^'5n< + 7— ^— Tg'^U + 7vT'^n< + ••• (5-3) 



n 



1 n — 1 (n — Ij^ (n — Ij'^ 



Clearly, the sum ends after a finite number of terms, since 5„< is a polynomial. 

Thus, the asymptotic solution (|5.ip is well-defined and uniquely determined as a function of 
the arbitrary constant b. To the lowest orders, we find 

If /3i = but /32 < we have the expansion 

1 ^Mlnt) 1 /33 , /3|-/32/?4 Int , 6 ^,0^.-2, ,A 
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If the first non- vanishing coefficient is /3„ then the expansion begins with t""*^/". 

Now we generalize this result to the case of more couplings. Consider a theory with s couplings 
9 — {9i} beta functions (ii{g)- To be specific, we assume 

°° 1 1 
= X] —fih-u9h ■ ■ ■ 9j„ = ^CijkQjgk + 0{g^), (5.6) 

n=2 

where the constants Cijk are the one-loop coefficients. We look for asymptotic solutions of the RG 
equations starting form the "Zimmermann trajectories" 

9^{t) - -J (5.7) 

in the limit t — )• cxd. Inserting (j5.7p into the RG equations 

m{t) = I3,{g{t)), (5.8) 



where Pi are given by (|5.6p . and keeping only the leading terms, we see that the constants at are 
determined by the quadratic equations 

1 

Oi = -Cijkdjak, (5.9) 

which have, in general, a discrete set of solutions. Normally the solutions just have to be real 
(if the couplings are parametrized to be real, as we assume), but in some cases further physical 
restrictions might apply. For example, stability (positive-definiteness of the action in Euclidean 
space) might require that some couplings be positive. We do not consider such restrictions here 
and assume that all real solutions are physical acceptable. It is straightforward to adapt our 
conclusions to specific situations. 

Around a Zimmermann trajectory, we continue the expansion as 

gi{t) ~ (ai + ^] , (5.10) 



t \ 

assuming that 7 is a positive number. Then 7 and bi are an eigenvalue and an eigenvector of the 
real matrix 

respectively. 

The matrix M is crucial for our discussion. The dimension ^af of the domain Paf of asymp- 
totic freedom is equal to the number of AA-positive eigenvalues 7 , including multiplicities. Observe 
that because of (j5.9p one eigenvalue is always equal to 1, with eigenvector Oj. If the Zimmermann 
trajectory exists, the dimension ^af is at least 1. If 7 = 1 is the unique positive eigenvalue, the 
form of the expansion is (|5.ip . 
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The most general solution reads 



Si(«) 




E 

7-n,7i>0 



6i,^.n(lnt) 



(5.11) 



Here 7 is a vector collecting the positive eigenvalues of the matrix A/", while n is a vector of 
non-negative integer entries. The condition n > means that n must not vanish identically. Two 
vectors n and n' such that 7 • n = 7 • n' are equivalent, and associated with a unique numerator 
&i,7-n- The sum is ordered for increasing values of 7 • n. Finally, the 6j.^.n(lnt)'s are polynomials 
of certain finite degrees in Int. 

Inserting (jS.lip into the RG equations (|5.8p and isolating the coefficients of the powers t"^"'"''"", 
we obtain equations for the polynomials bi^^.^- It is immediate to find that such equations have 
the form 



where (5j, 



7n< 



[(7 • n)(5ij 
is a sum of monomials 



b' 



i,-f-n 



-'i,7-n< ! 



(5.12) 



with 7 • iifc < 7 • n and J2k 7 ' nfc ^ 7 • n. 

Clearly, (^i,^.n< contains a finite number of terms. Now we want to show that equations (15. 12^ 
allow us to recursively determine the bi^-y.nS. 

For intermediate purposes, it is convenient to turn to the basis where the matrix A/ij has a 
real canonical Jordan form. Quantities in this basis are denoted with a tilde. Specifically, J\f is 
block-diagonal. Its first blocks are associated with the real eigenvalues A and have the form 

/ A \ 



1 

V 



(5.13) 



while the last blocks are associated with the complex eigenvalues fi and have the same form as 
(|5.13p . where however the A's are replaced by 2x2 blocks 

(Re/i Im^ \ 
— Im/i Re^u I ' 

the I's are replaced by 2x2 identity matrices and the O's are replaced by 2x2 matrices with 
vanishing entries. All matrices J\fij — (7 • n)dij are then in canonical Jordan form. 



Let Mij be such that J\f = W'^MM and h 
M to the left, we can rewrite it in the form 



«,7-n 



Mijbj^^.n- Multiplying equation (I5.12p by 



(7 • n)5ij - Mi^ 



b' 



«,7-n 



^i,-f-n< ) 



(5.14) 
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By induction, if we assume that the polynomials &j,^.n' with 7 • n' < 7 • n are known, we conclude 
that the 5j,^.n<'s are polynomials of certain finite degrees in Int. When the matrix Nij — (7 • n)6ij 
is invertible, the polynomials &j,^.n are uniquely determined. Calling lAij the inverse matrix of 
Nij — (7 • n)5jj, we have, similarly to (j5.3p . 



^i,7-n — \^'j,')-n< + ^j',7-nj ~ ~^ijSj,^-i[i< ^ ^ij^jk^'k,-y n< ~ ^iMjk^kl^l,'y-xi< + " " " (5.15) 

Again, the sum ends after a finite number of terms, since the (5i^'y.n<'s are polynomials. 

When J\fij — (7 • n)5jj is not invertible, one of its blocks has rriQ > zeros on the diagonal. 
Assume that this block is the one with i ^ i,j < i + rriQ and proceed as follows. The block- 
structure of A/" allows us to split equation (j5.14p into: a) the equation for the 6j.^.n's with i <r, b) 
the equation for the bi^^.^ with i ^ i, j < i + ttt-q; c) the equation for the 6j^^.n's with i ^ i + rriQ. 
Equations a) and c) are solved by formulas similar to (|5.15p . Equation b) has the form 



/ \ 

Ci ■ ■ • 

Cmo-i / 



^j,7-n ^4,7-11 — '^i,7-n<) 



where the QiS can be equal to or 1. The right-hand side is made of recursively known polynomials 
of some degrees di in Int. Then the &j,^.n's are polynomials of finite degrees greater than di, and 
each of them is uniquely determined up to an arbitrary additional constant. Therefore, in total 
we have mg arbitrary constants. 

Thus, equations (j5T2]) - ([5Tll) can be solved recursively to determine the polynomials ^2,7-n? 
and therefore the 6j.^.n's. The solutions (|5.1ip contain a number of arbitrary constants equal to 
the number of times the matrices Nij — (7 • n)5jj become degenerate, including multiplicities. 
This number is equal to the number of positive AA-eigenvalues, including multiplicities. Indeed, 
recalling that different n's with the same 7 • n are equivalent, each equation 7j =7-11 admits 
precisely one solution, and no other degeneracies are possible. 

The set of arbitrary constants contained in the asymptotic expansion (|5.1ip determines the 
domain of asymptotic freedom Paf- We conclude that the dimension of Paf is equal to the 
number of positive eigenvalues 7 of the matrix M, including multiplicities. 

In practice, we have to look for the Zimmermann trajectory around which the asymptotic 
expansion (|5.1ip has the maximal number of positive eigenvalues. In most cases the other Zim- 
mermann trajectories also play a role. They can determine the boundary of PaFj if it is two- 
dimensional, or its edges, if it more than two dimensional. 

Now we consider two examples in detail. 
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Dirac model First we consider the Dirac model (|4.6p in the large limit, focusing on the 32-54 
subsystem. The RG flow is given by the equations 



132 = 92 = K{gl + gl) + 0{g^), (3^ = g^ = 2^32^4 + 
where k = N/[18tt'^) > 0. We find the constants 

C222 = C244 = C424 = C442 = 2k, 

while all other entries Cjjfc vanish. The Zimmermann trajectories are given by (02,04) = (1,0)/k 
and (02)^14) = (1; Expanding around the trajectories with (02,04) = (1,0)/k we find 

that J\f is equal to the identity matrix, so 7 = 1 with multiplicity 2. This means that T>af has 
dimension 2. Two arbitrary constants appear at order 
Precisely, we find 

g2it) = -^^ + ^ia^nt + b2) + 0{t-HnH), g^{t) = 1(^4 Int + 64) + O (^"^ In' t) , (5.16) 

where 62,4 are arbitrary constants, while ,^2,4 are uniquely determined by the cubic terms of the 
beta functions. Since all other matrices J\fij — {'y-n)Sij are invertible the higher-order corrections are 
uniquely determined. Thus, the two-dimensional domain of asymptotic freedom is parametrized 
by the arbitrary constants 62,4- The other two Zimmermann trajectories, given by (02,04) = 
(1, ±1)7(2^), are the boundary of Paf- 

At one loop we can check our results solving the system explicitly. Call X± = 5^2 i 54- We 
have 

/3± = A± = kX^, 

wherefrom 

In the ultraviolet limit t — t- 00 the asymptotic behaviors are 

Here the logarithmic terms of (I5.16P are absent, because we are neglecting higher-order corrections 
to the beta functions (check also (j57 



Electroweak model Now we apply our method to the electroweak four fermion model (14. 4p 
with hi = bn = 1. We first restrict to the subspace {g2L^ A,ti). At one loop we have 

/32L = 77^2(352^ + 2A2), /3a = (2n - 5A) , /3„ = j^i^X^ + Qgl^ + 2Au). (5.18) 

(67r)"^ (67r)^ (67r)^ 
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The conditions j32L = f^x = Pu = ^ are solved by §21 = A = 0, while u remains arbitrary. Since 
a non-trivial fixed point cannot be trusted within our approximation, we need to project onto a 
suitable subspace of parameter space. Such a projection is automatic in our approach. 
The Zimmermann trajectories are 

ai a2 as 

52L ~ — -, A~ — -, u ~ — -, 

with 

(ai, 02, as) = 0, , Q, ^, 1^ , (0.046, -0.141, 0.149) . (5.19) 

The third trajectory is given by complicated irrational coefficients of which we just give the 
approximate numerical values. Two complex solutions to (j5.9p also exist, but must be discarded. 
Expanding around the first Zimmermann trajectory, we find 



1 





\ 





1 
3 





4 


4 
3 


-1/ 



(5.20) 



so the positive eigenvalues are 7 = 1/3 and 7=1. The Jordan canonical form of this matrix is 
diagonal: M =diag(l, 1/3, — 1). The domain of asymptotic freedom is two-dimensional and the 
arbitrary constants appear at orders and Using the one-loop truncated beta functions, 

the asymptotic expansion of the solution reads 

1 fin^ 1 1 5fifin^ ^n'^ 

(67r)-252L(t) = - - + ^ + ^(72a3 Int + 27 + b) - - ^(SGOa^ Int + 2907a3 + 56) 



3888a' 



t3 ' 

, 9 , , , a 9a^ 63a3 a , , , q , , , o 

(67r)-2A(t) = _ + _ + ___ ^(144a3 Int + 26 - 4950^ 



a2 / o 138267 o 45 \ 

- ^ (l620a=^lnt " -^a' + ybj , (5-21) 

(6.)-n(t) = - | + ^ + ^ + ^(72a3 Int + 6) - ^ (l44a3 Int + ^a^ + 26) 



3a^ / s, 70587 . 25 \ .In^t 
gOOa^ Int + —r—a^ + ^6 - 7776a^ 



t8/3 V 14 2 y t3 

up to 0{t~^lnt), where a and 6 are the arbitrary constants. 

The cubic corrections to the beta functions start with terms ~ l/t^, so they give extra contri- 
butions of the form ~ Cj/t^ to the solutions, where are uniquely determined functions of a and 
6. The corrections do not affect the terms proportional to (In t)/t^. Thus, the complete solution 
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has the form 



{Q^r^m = — + — + + 0(t-7/3 1^^)^ (5 22) 

2 1 'in^ 1 

(67r)-\(t) = -- + _ + -^ + -(72a3 Int + 6 + ^3) + ©(t"'/^ Int), 

where i = 1,2,3, are calculable numbers, depending on the cubic corrections to the beta 
functions. The other beta functions of (|4.4p give 

= -ICj + 1 + ^''^r'sn = + 1 + or"'). 

where ^4^5 are calculable numbers. 

We have thus found a two dimensional domain of asymptotic freedom. 

One-loop degeneracies In special cases, not frequent in physical problems, the one-loop coef- 
ficients Cijk can have degeneracies that make the expansions of some couplings start from powers 
^-1/" instead of 1/t, similarly to what happens in (|5.5p when /3i = for a single coupling. Some 
higher loop contributions can be as important as the one-loop ones, or even more important than 
the one-loop ones. Then the expansions of the beta functions in powers of the couplings have to 
be accordingly reordered. For example, consider the system 

91= 9i+ K^9igh 92 = ^9192, 

where k is a constant. Observe that the terms and K?gig2 are one- and two-loop, respec- 
tively. Nevertheless, they are equally important in the asymptotic expansion. The Zimmermann 
trajectories are 

to,,,) = (_l,„),(4±_L), ,5,23) 

to which we must add the line of fixed points (71 = 0. The procedure described above has 
to be applied using the trajectories (|5.23p . For example, expanding around the second pair of 
trajectories we find that the matrix N has eigenvalue 1 with degeneracy 2 and the expansions 
read 

= -U2 + 2^ In t--\+ 0(t-' In^ (), 

where a and h are the arbitrary constants. All higher-order terms are uniquely determined. The 
domain of asymptotic freedom is two-dimensional. 
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Another way to find X'af Here we give an alternative method that can be useful to determine 
the domain Paf of asymptotic freedom when the origin is an isolated fixed point (possibly after 
a suitable projection in parameter space). We define the radius p in parameter space as 



N 



\ i=i 



and the radial velocity v as 



dt p^ 



Let V denote a domain in parameter space and V its closure. A theory is asymptotically free 
in Paf if the origin belongs to Paf (but not Paf) and every trajectory passing through Paf 
remains in Paf and flows to the origin in the ultraviolet limit. The trajectories that satisfy v <0 
asymptotically for t — > oo in a neighborhood of the origin belong to Paf- The Zimmermann 

trajectories, in particular, belong to Pap- 

When the origin (7 = is an isolated fixed point Paf can be also characterized as follows: 

i) find the domain T> around the origin where v < 0; 

ii) call d^D the boundary of T) minus the origin, and consider the trajectories crossing it; 
iv) if all such trajectories enter P, then Paf = if not, Paf is minus the trajectories 

leaving T> through the boundary d^D. 

Indeed, the remaining trajectories cannot leave P, so they must flow to the origin, because it 
is the unique fixed point. 

The condition i; < is necessary, but not sufficient, because some trajectories intersecting T> 
can cross its boundary and run away, instead of flowing to the origin. The good feature of T> is 
that it can be easily determined, but Paf is only a subset of T). 

Now, observe that T> depends on the parametrization of the couplings, while Paf of course 
does not. Call hi{g) a reparametrization of the couplings and P/j the domain where the velocity 



is negative in a neighborhood of the origin. An efficient way to estimate Pap (and in most cases 
determine it) is to take the intersection of the P?j's, for all rcparamctrizations h. 

We illustrate this method in the (72-54 subsystem of the Dirac model in the large N limit, 
using one-loop truncated beta functions. Observe that the origin is the unique fixed point. We 
have 




pv = g2P2 + 54/34 = 052(52 + 354)- 
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The domain 2? is just g2 < 0. Consider now the one-parameter family of reparametrizations 

^2=52+054, ^4 = 052 +54- 

We find 



PhVh = 




where r = 52/54 • Now we study the condition Vh < 0, knowing that 52 must be negative. Varying 
a to obtain the best result, we find 

\r\ < 1, i-e. I52I < |c/4|, 

which together with 52 < gives our best estimate of the domain T>af- 

We can check this estimate using the explicit one- loop solution (j5.17p . We see that inside V 
only the trajectories with 152! ^ I94I flow to the origin. All others cross the boundary doD (that 
is the line 52 =0), enter the region with v > and run away. Thus, Paf is given by 52 < 0, 
I52I ^ \94\- We conclude that our method gives an accurate estimate of PaFj since it misses only 
its boundary, namely the trajectories with \g2\ = \g4,\- 



6 Zimmermann trajectories and hidden symmetries 

In the previous section we have seen that the Zimmermann trajectories play an important role in 
the study of the domain of asymptotic freedom. Some such trajectories (see for example (|5.19p ) 
involve only rational coefficients, others very complicated irrational numbers. Normally, rational 
trajectories appear when the theory has more symmetries. In this section we point out that 
the existence of rational Zimmermann trajectories appears to be a general feature of high-energy 
Lorentz violating four fermion models. 

First, we briefly recall Zimmermann's "reductions of couplings" OE]. Assume that a theory 
has couplings A/, / = 1, . . . A^. Zimmermann's idea is to parametrize the couplings in terms of 
a smaller set of independent parameters aj, j = 1, . . . M < A^. Write A/ = Xj{aj). Consistence 
with the renormalization group demands 

The most interesting case is A^ = 2, M = 1. Normally, if the one- loop beta functions are quadratic 
in the couplings, as in our case, Zimmermann's equations admit two power-series solutions of the 
form 

00 

\ = ca + a dno"', 

n=l 
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where c and are calculable, generically irrational, numbers. If the solution exists at one loop 
(namely, if c is real), then it exists to all orders. The most general solution to Zimmermann's 
equations is not analytic, but has the form 



A' = A+ d, 



m,n=l 



where ^ is typically irrational, dn is arbitrary and the other coefficients dmn ^ire uniquely deter- 
mined once dii is given. 

In special situations, such as when the "reduced" model has additional symmetries, the power- 
series solution A shrinks to the monomial ca, with a rational coefficient c, and A' becomes analytic. 

For example, the (Lorentz invariant) model of a spinor ip and a pseudoscalar field A interacting 
with the Lagrangian 



admits reductions 



A' = ^(1 ± ^145)5^ + di/ + • • • + + • • • 

On the other hand, the massless model with interaction 

Ci = g^{A + ij5B)i, - ^{A^ + BY 

admits the rational reduction 

A=5^ (6.1) 

which reveals the existence of a symmetry. Indeed, when the couplings are related by formula 
(|6.ip we have the supersymmetric Wess-Zumino model. 

We now analyze Zimmermann's trajectories at one loop in some of our Lorentz violating four 
fermion models. Consider again the electroweak model. We look for RG trajectories where all 
couplings are proportional to one another: gn = ax, g2L = bx, gji = cx, giR = dx and A = ex. 
The constants a, 6, c, d and e can be worked out matching the beta functions (j4.4p . The absolute 
value of 6/, (or 5^) can be set to 1 rescaling the space coordinates. On the other hand, the ratio 
bi/bn is free and does not run at one loop. 

We choose {bil = \bji\ = 1 and consider the cases s = ±1. We find only three real solutions. 
One of them is just gn = gn = giR = A = 0, with only g2L non-vanishing. Of the other two 
solutions, only one has rational coefficients, and reads 

giL = -, 52L = A, gR = -\ gLR= A. (6.2) 

5 5 lU 

If we choose ^ we generically do not find rational trajectories. 
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The existence of a trajectory with rational coefficients is unexpected and offers evidence that 
the reduced model might have hidden symmetries. Its Lagrangian reads 

_ 9 |r]iri^2 9 (T]iTj\lT]j Ti\ '^^ , \ ^ .2 (5 + Ss)^ 



\^kl^ > IKY ' 5Ai'^^' 10A2 

where g = —A > 0, and its one-loop beta function is 

5g2 



(67r)2 • 



If the hidden symmetry is simple, we expect that the relations (|6.2p are preserved by higher- 
loop corrections. However, this is not a necessary requirement for a hidden symmetry. 

The rational RG trajectory exists also in the U{Nl) x U{N{i) model for = \hpt\ = 1. We 
find 

Nl + Nr-2 2Nl + 2Nr + 3s-1 

5iL-5ii?--2^^^2Ar^_-^9, 92L - g2R - 92LR - -g, giLR - - 2{2Nl + 2Nr-1) ^' 

the beta function being 

^g = -^{2NL+2NR-l). 

Curiously, the model of Dirac fermions in the large N limit admits only rational trajectories. 
We have already studied the 92-54 subset in the previous section. In the 55,6-/1,2 subset we find 
15 rational trajectories: 

55 = (-2, ±1, 1, -1, -5, -2, -2, -7/5)56, 

/i = (0, 0, ±1/2, ±1/2, ±3/2, ±1, ±1/2, ±1/10)56 , 

/2 = (0, 0, 0, ±1, ±3, ±3/2, ±3/2, 0)56 • 

One such trajectory (55 = ge = 2/i, /2 = 0) gives the Lorentz violating "Gross-Neveu" model 
[lOj . whose interaction reads 

which is renormalizable by the same argument used in (|4.2p . In some sense, this is an example of 
"hidden symmetry" associated with the rational trajectory. 
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7 Conclusions 



If Lorentz symmetry is violated at high energies, then the Standard Model admits a ultraviolet 
completion that is renormalizable despite it contains four fermion vertices at the fundamental 
level. In this scenario, four fermion models play a key role, because all other interactions, being 
super-renormalizable, disappear at energies much higher than the scale of Lorentz violation. In 
this paper we have studied the one-loop renormalization of CPT-invariant Lorentz violating four 
fermion models and their RG flows. 

We have first considered the most general case, working out formulas for the beta functions, 
and then analyzed particular models in detail. We have formulated a method to determine the 
domain of asymptotic freedom expanding the running couplings around the free fixed point. We 
emphasize that if the four fermion model (II. 3p is asymptotically free then the entire Lorentz 
violating Standard Model is. 

Moreover, we found that the RG flow admits a number of special "rational" trajectories that, 
in the spirit of Zimmermann's reduction of couplings, might hide some new symmetries. 

References 

[1] V.A. Kostelecky and N. Russell, Data tables for Lorentz and CTP violation, V.A. Kostelecky, 
Ed., Proceedings of the Fourth Meeting on CPT and Lorentz Symmetry, World Scientific, 
Singapore, 2008, p. 308 and larXiv:0801.0287l [heo-ph]. 

[2] D. Anselmi and M. Halat, Renormalization of Lorentz violating theories, Phys. Rev. D 76 
(2007) 125011 and arXiv:0707.2480 [hep-th]. 

[3] D. Anselmi, Weighted power counting, neutrino masses and Lorentz violating extensions of 
the Standard Model, Phys. Rev. D 79 (2009) 025017 and >arXiv:0808.3475, [hep-ph] . 

[4] D. Anselmi, Standard Model Without Elementary Scalars And High Energy Lorentz Viola- 
tion, Eur. Phys. J. C 65 (2010) 523 and iarXiv:0904. 18491 [hep-ph]. 

[5] W. Zimmermann, Reduction in the number of coupling parameters, Commun. Math. Phys. 
97 (1985) 211. 

[6] See for example, A. Strumia and F. Vissani, Neutrino masses and mixings and..., 
arXiv:hep-ph/0606054] [hep-ph] . 

[7] W. Zimmermann, Scheme independence of the reduction principle and asymptotic freedom 
in several couplings, Commun. Math. Phys. 219 (2001) 221. 



23 



[8] R. Oehme, Reduction and reparametrization of quantum field theories, Prog. Theor. Phys. 
Suppl. 86 (1986) 215. 

[9] R. Oehme and W. Zimmcrmann, Relation between effective couplings for asymptotically free 
models, Commun. Math. Phys. 97 (1985) 569. 

[10] D.J. Gross and A. Neveu, Dynamical symmetry breaking in asymptotically free field theories, 
Phys. Rev. D 10 (1974) 3235. 



24 



